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Using a 478 pb _1 data sample collected with the BESIII detector operating at the Beijing Electron Positron 
Collider storage ring at a center-of-mass energy of yfs = 4.009 GeV, the production of e + e~ — > rjj /ip 
is observed for the first time with a statistical significance of greater than lOcr. The Born cross section is 
measured to be (32.1 ± 2.8 ± 1.3) pb, where the first error is statistical and the second systematic. Assuming 
the r/J/tp signal is from a hadronic transition of the ?/>(4040), the fractional transition rate is determined to 
be S(V>(4040) -> ■qj/ip) = (5.2 ± 0.5 ± 0.2 ± 0.5) x 10~ 3 , where the first, second, and third errors are 
statistical, systematic, and the uncertainty from the i/>(4040) resonant parameters, respectively. The production 
of e + e _ — > 7T° J/ip is searched for, but no significant signal is observed, and S(t/j(4040) — > Tv°J/ip) < 
2.8 x 10~ 4 is obtained at the 90% confidence level. 
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PACS numbers: 13.25.Gv, 13.40.Hq, 14.40.Pq 

The properties of excited J PC = 1~~ charmonium states 
above the DD production threshold is of great interest but not 
well understood, even decades after their first observation HI. 
The current experimentally well established structures in the 
hadronic cross section are the 7/>(3770), '0(4040), "0(4160), 
and 0(4415) resonances 0]. Unlike the low-lying vector cc 
states J/ip and -0(3686), all of these states couple to open- 
charm final states with large partial widths, and disfavor hid- 
den charm decays. 

Recently, new vector charmonium-like states, the y(4260), 
the F(4360) and the F(4660) have been discovered via their 
decays into exclusive n + TT~J/ip and ir + Tr~ 0(3686) final 
states Jit]. The common properties of these states are rela- 
tively narrow widths and strong couplings to hidden-charm 
final states. These F-states cannot be assigned to any of the 
conventional cc family states J3l in any natural way 
and suggest the existence of a non-conventional meson spec- 
troscopy 

Hadronic transitions play an important role in understand- 
ing the nature of conventional heavy quarkonium. An ex- 
cess of 77 over 7T + 7r _ hidden-bottom transition rates of the 
T(45) Jg] has been explained as an admixture of a four- 
quark state in the T(45) wave function |7]. A similar picture 
might be expected in the charm sector but, as of yet, there is 
no experimental data available for r/ transitions in the high- 
mass charmonium and charmoniumlike states, except for evi- 
dence of 0(3770) rjJ/ip (3.5a) H and 0(4160) rjj/tp 
(4.0a) fl. Moreover, there are predictions of many new states 
in various models trying to explain the conventional and un- 
conventional states observed in this mass region Jj§]. 

In this Letter, we report cross section measurements for 
e + e~ — > r]J/i[> and ir°J/ip at the center-of-mass energy 
\/s = (4.009 ± 0.001) GeV. The analysis is performed with 
a 478 pb _1 data sample collected with the BESIII detector 
located at the BEPCII storage ring iflQll . The integrated lumi- 
nosity of this data sample was measured using Bhabha events, 
with an estimated uncertainty of 1 . 1 %. In order to control sys- 
tematic errors, an accompanying data sample of about seven 
million 0(3686) events was accumulated under the same ex- 
perimental conditions. In the analysis, the J/ip is recon- 
structed through its decays into lepton pairs (e + e~~ and 
while 77/tt is reconstructed in the 77 final state. 

The GEANT4-based Monte Carlo (MC) simulation soft- 
ware, which includes the geometric description and the de- 
tector response, is used to optimize the event selection crite- 
ria, determine the detection efficiency, and estimate the back- 
grounds. Signal e + e~ — > r/J/ijj and Tr°J/ipMC samples con- 
taining 20,000 events for each channel are generated. Initial 
state radiation (ISR) is simulated with KKMC lUHl . assuming 
ijj/ip and 7T° J/0 are produced via 0(4040) decays, and the 
0(4040) is described by a Breit-Wigner (BW) function with a 
constant width. The maximum energies of the ISR photons are 
347 MeV and 700 MeV, corresponding to 77 J/0 and ir°J/ip 



production thresholds, respectively. For backgrounds studies, 
MC samples equivalent to 1 fb _1 integrated luminosity are 
generated: inclusive 0(4040) decays, ISR production of low- 
mass vector charmonium states, and QED events. The known 
decay modes of the charmonium states are generated with 
EVTGEN lfl2ll with branching fractions set to their world av- 
erage values lull an d the remai ning events are generated with 

LUNDCHARM |Qj] or PYTHIA 11411 ■ 

Charged tracks are reconstructed in the main drift chamber, 
and the number of good charged tracks is required to be two 
with zero net charge. For each track, the polar angle must sat- 
isfy I cos#| < 0.93, and the point of closest approach to the 
e + e~ interaction point must be within ±10 cm in the beam 
direction and within ±1 cm in the plane perpendicular to the 
beam direction. A charged track with deposited energy in the 
electromagnetic calorimeter less than 0.4 GeV is identified as 
a n candidate while that with a deposited energy over mo- 
mentum (E/p) ratio larger than 0.8 is identified as an electron 
candidate. Both of the two charged tracks are required to be 
either identified as muons or as electrons. 

Showers identified as photon candidates must satisfy fidu- 
cial and shower-quality requirements. The minimum energy 
is 25 MeV for electromagnetic calorimeter barrel showers 
(| cos0| < 0.8) and 50 MeV for end-cap showers (0.86 < 
I cos#| < 0.92). To eliminate showers produced by charged 
particles, a photon must be separated by at least 20 degrees 
from any charged track. Final state radiation (FSR) and 
bremsstrahlung energy loss of leptons are corrected by adding 
the momentum of photons detected within a 5 degree cone 
around the lepton momentum direction. The number of good 
photon candidates is required to be two (the efficiency is over 
95%), and the recoil mass of the two photons M reco ii(77) = 
^{Pcm - PI- P2) 2 e [2.9, 3.4] GeV/c 2 is required to se- 
lect good J/0 candidates. Here Pcm is the four-momentum 
of the initial states, and PI, P2 are the four-momenta of the 
two photons. 

The lepton pair and the two photons are subject to a four- 
constraint (4C) kinematic fit to improve the momentum reso- 
lution and reduce the background. The chi-square (x 2 ) of the 
kinematic fit is required to be less than 40. In order to reject ra- 
diative Bhabha and radiative dimuon (-fe + e~ /jfi + fi~) back- 
grounds associated with an energetic radiative photon (jh) 
and a low energy fake photon, the invariant mass M (^H^ + i~ ) 
is determined from a three-constraint (3C) kinematic fit in 
which the energy of the low energy photon is allowed to 
float. Since the fake photon does not contribute in the 3C- 
fit, the M(7^f£ + ^ _ ) mass distribution is not distorted by the 
photon energy threshold cutoff, and backgrounds are clearly 
separated from signal. The requirement M(jh£ + £~) < 
3.93 GeV/c 2 removes over 50% of radiative Bhabha and ra- 
diative dimuon background events with an efficiency greater 
than 99% for 77 J/ip and 89% for n°J/tp. 

After imposing all of these selection criteria, the invariant 
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mass distribution of lepton pairs is shown in Fig. Q] A clear 
J/i[> signal is observed in the mode while indications 

of a peak around 3.1 GeV/c 2 also exist in the e + e" mode. 
The remaining dominant backgrounds are surviving radiative 
dimuon events in /J + /i~ and radiative Bhabha events in e + e~ ; 
these contribute fiat components in the M{£ + £~) distributions 
with no associated peaks in the M^j) invariant mass dis- 
tribution. The high background level in the e + e~ mode is 
due to the huge background from the Bhabha process. Other 
possible background sources include e + e~ — > ir a ir Q J/tp, 
7r+7r-7r°/7r+7r-?7, and 7x c j(lP)/7Xcj(2P). The 7r°7r° J/ip 
background is estimated by MC simulation to be at the 4.5 pb 
level and, thus, negligibly small |@]. Potential jXcjO-P) and 
7Xcj(2-P) radiative transition backgrounds are estimated us- 
ing the selected data sample; no significant signal is found 
for either XcjO-P) or x c j(2P) in M(jJ/ip) mass distribu- 
tion. The 7T + 7r~7r° and ir + 7r~rj backgrounds are estimated 
using J/ip sideband events. The ISR-produced vector char- 
monium backgrounds, including JisrJ/iP, 7isr' ( /'(3686) and 
7isrV'(3770), are estimated by means of an inclusive MC 
sample and only 3.3 events in the /-t + /i~ mode and 3.1 events 
in the e + e~ mode are found (normalized to data luminosity). 
As they would peak at neither the 77 nor the ir° signal region, 
they are neglected in the analysis. 

The resolution of the invariant mass of the lepton pairs is 
determined to be 14 MeV/c 2 by MC simulation, and is in 
good agreement with events in the -0(3686) data sample. The 
mass window of the J/ip signal is defined as 3.075 GeV/c 2 

< M (£+£-) < 3.125 GeV/c 2 , and the sidebands are defined 
as 2.95 GeV/c 2 < M(£+£~) < 3.05 GeV/c 2 or 3.15 GeV/c 2 

< M{£ + £~) < 3.25 GeV/c 2 , which is four times as wide 
as the signal region. Figure |2] shows the M(77) invari- 
ant mass distributions for events in the J/ip —> and 
J/ip —> e + e~ signal regions. A significant 77 signal is ob- 
served in both modes. In the M(77) distribution for J/ip 
mass-sideband events, there are backgrounds that peak in 
the 7T° signal region in J/ip — > [i + [i~ that originate from 
e + e _ — > 7T + 7T~7r . In order to suppress e + e~ — > 7r + 7r~7r° 
backgrounds, at least one charged track is required to have a 
muon counter hit depth larger than 30 cm for the TT°J/tp signal 
search. The efficiency for this requirement is 87.9% for sig- 
nal while about 74% e + e~ — > 7r + 7r~7r° background events 
are rejected. Figure [3] shows the M (77) invariant mass distri- 
bution below 0.3 GeV/c 2 for J/ip -4 /i + /J~. No significant 
7T° signal is observed. We do not analyze ir°J/ip production 
in J/ip — > e + e~ due to the huge background from Bhabha 
events. The final selection efficiencies are 38.0% in /i + /J~ 
and 26.9% in e + e _ for r) J/i/j, and 31.1% in for ir° J/ip, 
according to MC simulation. 

The 71/(77) invariant mass distributions are fitted using 
an unbinned maximum likelihood method for M(7*y) < 
0.9 GeV/c 2 in both modes. The probability density function 
(pdf) for the 77/71-° signal in J/tp — > is taken from MC 

simulation, while in J/tp — > e + e~, only the 77 pdf from MC 
simulation is used. To account for resolution differences be- 
tween data and the MC simulation, three Gaussian functions 



are convolved with the r\ and the 7T° signal pdfs. For the 
77 signal, the standard deviation of these Gaussians are free 
while for tt° signal, it is fixed to (2.4 ± 0.9) MeV/c 2 , which 
is determined from a -0(3686) — > ir°J/ip control sample. 
Background shapes are described by a third-orderpolynomial. 
Figure [2] shows the fit results for the 77 signal and the back- 
ground contributions for J/ip — > and J/ip — > e + e~. 
The fits yield (77) = 111.4 ± 11.0, and iV^ e _ (77) = 
61.4 ± 10.5. The standard deviation of the smearing Gaus- 
sian convolved with the 77 signal is (3.7 ± 1.0) MeV/c 2 in 
and (3.7 ± 1.9) MeV/c 2 in e + e _ . Good agreement 
is observed between the two modes, and these values are 
consistent with values from the t/;(3686) — > r]J/ip control 
sample (3.4 ± 0.6 MeV/c 2 in and 4.6 ± 0.6 MeV/c 2 
in e + e~). The goodness of fit is estimated by using a \ 2 
test method with the data distributions regrouped to ensure 
that each bin contains more than 10 events. The test gives 
X 2 /W./=14.1/14=l.l fotfi+fj,- and x 2 /W/=42.9/43=1.0 
for e + e~. Figure [3] shows the fit result for the 7r° signal and 
the background contribution for J/ip — > /x + . Since the 7r° 
signal is not significant, we determine an upper limit for the 
7T° signal yield of AT u p(7t°) < 11.7 at the 90% confidence 
level. The e + e~ — > ir + ir~ir° backgrounds are estimated 
by fitting the M (77) distribution of the J/tp mass sideband 
events. The signal pdf for the ir° is a Gaussian function and 
that for the background is a third-order polynomial. The fit 
yields N + s _ (ir°) = 2.8 ± 1.1 after normalization. The sta- 
tistical significances of the 77 and ir° signals are examined by 
means of the difference in log-likelihood value with or with- 
out signal in the fit and the change of the number of degrees of 
freedom (Andf). For the 77 signal, the statistical significance 
is larger than 10ct while that for the ir° signal is only 1.1c 

The Born-order cross section is determined from the rela- 
tion 

B _ N Rt - N hk s 
a ~ £ int (l + <J)eB' ( } 

where N Rt and A^ bkg are the number of signal events from the 
fit and the number of background events, respectively; C- int is 
integrated luminosity; e is selection efficiency; B is branch- 
ing fraction of intermediate states decay; and (1 + 5) is the 
radiative correction factor, which is 0.757 according to QED 
calculation [15]. 

For the e + e~ — > r]J /ip cross section, we obtain a B = 
34.8 ± 3.5 pb for the fj+n~ mode, and a B = 27.1 ±4.7 pb for 
the e + e~ mode. Since the results from the two modes agree 
with each other, we quote a combined cross section result: 

a B {e + e- -^r]J/ip) = 32.1 ± 2.8 pb. (2) 

Here the errors are statistical only. 

Systematic errors mainly come from the luminosity mea- 
surement, detection efficiency, background estimation and 
branching fractions of intermediate states decays. All the con- 
tributions are summarized in Table U 
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FIG. 1 : (Left panel) M([i + fi ) and (right panel) M(e + e ) invariant mass distributions. Dots with error bars are data and the open histogram 
in the left panel shows inclusive-MC-estimated background events. 




FIG. 2: Distributions of M(pfy) between 0.2 GeV/c 2 and 0.9 GeV/c 2 for J/rfi — > (left panel) and for J/ip — > e + e~ (right panel). Dots 

with error bars are data in J/ip mass signal region, and the green shaded histograms are from normalized J/xp mass sidebands. The curves 
show the total fit and the background term. 




FIG. 3: Distribution of M(77) below 0.3 GeV/c 2 for J /ip — > . Dots with error bars are data in J/ip mass signal region, and the green 

shaded histogram is from normalized J/tfi mass sideband. The curves show the total fit and the background term. 
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TABLE I: Summary of the systematic errors (%) in the cross section 
measurement. 



Source 


+ - 


rje e 


u + - 

TV )J, fJ, 


Luminosity 


1.1 


1.1 


1.1 


Tracking 


2 


" 


2 


Photon detection 


2 


2 


2 


Lepton resolution 


1.6 


2.4 


1.6 


Kinematic fit 


1.9 


1.9 


1.9 


Background shape 


1.5 


3.0 


9.4 


Fit function 






3.9 


V>(4040) parameters 


2.0 


3.3 


4.0 


Branching fractions 


1.2 


1.2 


1.0 


Others 


1.0 


1.0 


1.0 


Total 


5.0 


6.1 


11.8 



The uncertainty from luminosity measurement is estimated 
to be 1.1% using Bhabha events. The muon tracking effi- 
ciency is estimated to be 1 % for each track. Since the luminos- 
ity is measured using Bhabha events, the tracking efficiency 
of electron pairs cancels. The photon detection efficiency is 
also estimated to be 1% for each photon. The uncertainties 
associated with the lepton pair invariant mass resolutions and 
the kinematic fits are estimated using the -0(3686) — > 77J/0 
control sample. It is obtained from the -0(3686) data sample 
by imposing the selection criteria described above, and requir- 
ing M(-fn J/ip) < 3.49 GeV/c 2 to reject \ c \ and \c2 events. 
This gives a low-background 0(3686) — > 77 J / events with a 
purity of 98.5%. The efficiency difference between data and 
MC simulation for the J/ invariant mass window is 1 .6% in 
the /Lt + /i~ mode and 2.4% in the e + e~ mode. They are taken 
as systematic errors due to lepton-pair invariant mass resolu- 
tion. For the kinematic fit, the efficiency difference between 
data and MC simulation is 1.9% in both modes. 

Uncertainties due to the choice of background shape are es- 
timated by varying the background function from a 3rd-order 
polynomial to a 2nd-order and a 4th-order polynomial in the 
fit, and these changes yield a 1.5% difference in fi + fi~ and 
a 3.0% difference in e + e~ in the number of r\ signal events. 
The e + e~ — > 7r + 7r~7r° backgrounds subtraction gives a 9.4% 
difference in /i + /J~ in the number of 7r° signal events. The 
uncertainty due to the fit function is estimated by changing 
the smearing Gaussian parameter by one standard deviation 
in the 7r° signal pdf, which gives 3.9% difference in the num- 
ber of 7T° signal events. Uncertainties in the 0(4040) reso- 
nance parameters and possible distortions of the 0(4040) line 
shape due to interference effects with the nearby 0(4160) res- 
onance introduce uncertainties in the radiative correction fac- 
tor and the efficiency. Changing the Breit-Wigner parame- 
ters (mass and width) by one standard deviation according to 
PDG values H], or using a coherent shape with the 0(4160) 
resonance 1 1 611 result in variations in (1 + 5) x e of 2.0% in 
fi + \i~ and 3.3% in e + e~ for the 7/J/0 measurement, and, 
4.0% in for ir° J/ip measurement. The PDG uncer- 

tainty in B(J/ip (.+£-) is 1% and B(i] 77) is 0.5% Q. 



Other sources of systematic error, including fake photon sim- 
ulation and the final-state radiation simulation, are estimated 
to be 1.0% in total. 

Assuming all the sources are independent, the total system- 
atic errors on the 77J/0 cross section measurement is deter- 
mined to be 5.0% for /J + /i~ and 6.1% for e + e~. Considering 
the common and uncommon errors for these two modes, the 
combined systematic error on the 7/J/0 cross section mea- 
surement is 4.0%. The total systematic error is 11.8% in 
for the 7T° J/0 cross section measurement by summing 
up all the errors in quadrature. 

Since the significance of the 7r°J/0 signal is low, an up- 
per limit on the 7r°J/0 production cross section is set at 
a B (e + er ir° J/0) < 1.6 pb at the 90% confidence level, 
where e + e~ — > ir + TT~ir° backgrounds have been subtracted 
and the efficiency is lowered by a factor of (1 — a sys ). 

If we assume the observed 77 J/0 and 7r°J/0 are com- 
pletely from 0(4040) decays and use the total cross section 
of 0(4040) at y/s = 4.009 GeV [(6.2 ± 0.6) nb] calculated 
with the PDG resonance parameters J2l as input, we deter- 
mine the fractional transition rate #(0(4040) — > 77 J/0) = 
(5.2 ± 0.5 ± 0.2 ± 0.5) x 10~ 3 , where the first, second, and 
third errors are statistical, systematic, and uncertainty from 
0(4040) resonant parameters, respectively. In addition, we 
obtain an upper limit on S(0 (4040) tt°J/0) < 2.8 x 10~ 4 
at the 90% confidence level. 

In summary, we observe for the first time e + e~ — > 77 J/0 
production at = 4.009 GeV with a statistical significance 
greater than 10er. The Born cross section is measured to be 
(32.1 ± 2.8 ± 1.3) pb, where the first error is statistical and 
second systematic. We do not observe a significant e + e~ — > 
7T° J/0 signal, and the Born cross section is found to be less 
than 1.6 pb at the 90% confidence level. These measurements 
do not contradict the upper limits set by CLEO experiment |0] . 
The 77 J 1 cross section measurement is within the range of 
the theoretical calculation and the 7r° J/0 upper limit does not 
exclude the prediction IU7I1 . A transition rate of 5 x 10 ~ 3 level 
is measured for 0(4040) — > r]J/ip, corresponding to a partial 
decay width at the 400 keV level, which is much larger than 
that for 0(3770) — > 77 J/0 JH and is more than two times of 
that for 0(4040) tt+tt" J/0 U. 
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